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Abstract 

The reduced (in the angular coordinate (p) wave equation and Klein-Gordon 
equation are considered on a Kerr background and in the framework of C°-semigroup 
theory. Each equation is shown to have a well-posed initial value problem, i.e., to 
have a unique solution depending continuously on the data. Further, it is shown 
that the spectrum of the semigroup's generator coincides with the spectrum of 
an operator polynomial whose coefficients can be read off from the equation. In 
this way the problem of deciding stability is reduced to a spectral problem and a 
mathematical basis is provided for mode considerations. For the wave equation it is 
shown that the resolvent of the semigroup's generator and the corresponding Green's 
functions can be computed using spheroidal functions. It is to be expected that, 
analogous to the case of a Schwarzschild background, the quasinormal frequencies of 
the Kerr black hole appear as resonances, i.e., poles of the analytic continuation of 
this resolvent. Finally, stability of the background with respect to reduced massive 
perturbations is proven for large enough masses. 



1 Introduction 



The stability of the Schwarzschild black hole was demonstrated by Kay and Wald [14] 
who showed the boundedness of all solutions of the wave equation corresponding to C°° 
data of compact support. Their proof rests on the positivity of the conserved energy. 



The problem is more subtle for Kerr space time. A conserved energy exists, but the 
energy density is negative inside the ergosphere. Hence the total energy could be finite 
while the field still might grow exponentially in parts of the spacetime. Papers by Press 
and Teukolsky |[22| , Hartle and Wilkins |§, and Stewart make the absence of expo- 



nentially growing normal modes very plausible. Whiting [31] has proven that there are no 
such modes, and in his proof he showed that normal modes grow at most linearly in time. 
Recent numerical evolution calculations ITBL O for slowly and fast rotating Kerr black 



holes show no sign of exponential growth. In the case of massive scalar perturbations of 
Kerr results of Damour,Deruelle and Ruffini ||, Zouros and Eardley |34|], and Detweiler 
point to the existence of unstable modes. These modes are very slowly growing with 
growth times similar to the age of universe. This fact complicates the numerical detection 
of such modes. 



Here we consider the reduced (in the angular coordinate ip) wave equation and Klein- 
Gordon equation on a Kerr background and in the framework of C°-semigroup theory. 
For this the mathematical framework from || is used. For each equation it is shown that 
the initial value problem is well-posed, i.e., has a unique solution which depends contin- 
uously on the data. Further, it is shown that the spectrum of the semigroup's generator 
coincides with the spectrum of an operator polynomial whose coefficients can be read off 
from the equation. In this way the problem of deciding stability is reduced to a spectral 
problem. For the wave equation it is shown that the resolvent of the semigroup's genera- 
tor and the corresponding Green's functions can be computed using spheroidal functions. 
It is to be expected that, analogous to the case of a Schwarzschild background, the quasi- 
normal frequencies of the Kerr black hole appear as poles of the analytic continuation of 
this resolvent. Finally, the stability of the background with respect to reduced massive 
perturbations is proven for large enough masses. This is done by applying an abstract 
stability criterium from ||. 

The Kerr metric in Boyer-Lindquist coordinates t, r, 9, p is given by 

/ 2Mr\ )2 4Marsin 2 0, , £,„ „ t/l2 
g= ( 1 — J dt 2 + dtdp - —dr 2 - T,d9 2 

( 2 2 2Ma 2 rsin 2 \ . „ , 2 
- r 2 + a 2 H )sm9d<p , (1) 



where M is the mass, a e [0, M\ is the rotational parameter, 

A := r 2 — 2Mr + a 2 , £ := r 2 + a 2 cos 2 9 . (2) 

The coordinates are constrained by— oo<t<+oo,r + <r<+oo,— n < p < n and 
< 6 < 7i where 

r+ := M + VM 2 -a 2 . (3) 

As a little reminder on the Kerr geometry we give the following basic facts relevant for 
the discussion of the wave equation. The coordinate vector field d/dr becomes singular at 
r = r+. This value of the radial coordinate marks the event horizon for Kerr spacetime. 
The coordinate vectorfield d/dt is null on the ergosphere 



r = M + VM 2 - a 2 cos 2 9 , (4) 

is spacelike inside and timelike outside. So t is not a time coordinate inside the ergosphere 
and therefore one might think that Boyer-Lindquist coordinates are unsuitable for a sta- 
bility discussion. It turns out that this is not the case for the methods (from semigroup 
theory) of this paper. Finally, the Kerr metric is globally hyperbolic outside the horizon 
and hence the Cauchy problem for the scalar wave equation is well posed for data on 
any Cauchy surface. This result is not used in this paper. Existence, uniqueness and 
continuous dependence of the solutions on the initial data is proved here, too. 



The reduced wave equation governing solutions of the form ip(t, r, 9, 0) = exp(irmp)u(t, r, 9), 



where m runs through all integers, is given by 
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A first inspection shows that 
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(6) 



on Q := (r + , oo) x (0, n). Hence the coefficient multiplying idu/dt is real-valued, strictly 
positive and bounded. Moreover the coeffcients multiplying the derivatives in r and 9 are 
real- valued and vanish at the horizon. As consequence (j5|) is singular at all points on the 
boundary offl. 



The structure of this paper is as follows. Section 2 contains the used conventions. Sec- 
tion 3 gives an initial value formulation for (H). There the equation is interpreted as an 
abstract equation 

((A + C)u + iBv) ^ 

for a differentiate function *(-u,t>) assuming values in an appropriate Hilbert space Y 
and in particular in the domain of a linear operator G. Here ' denote a Hilbert space 
derivative. The linear operators A + C and B will be read off from (^). B is the maximal 
multiplication operator in X given by the function multiplying idu/dt. The auxiliar 
operator C is a suitable negative multiple of the identity operator on X. The definition 
of A is more involved. A preliminary form A + C of A + C is given by the differential 
operator enclosed in square brackets. It has as a domain all complex-valued functions on 
Q which are twice continuously differentiate and have a compact support in Q. Moreover 
X is chosen such that Aq + C is symmetric. It will be obvious that this operator is 
in addition semibounded (from below). In the next step A + C will be defined as the 
Friedrichs extension of A + C. Using this existence, uniqueness, continuous dependence 
of the solutions on the initial data (|7|) follow from abstract theorems derived in B. In 
Section 4 the domain of A will be investigated further. This is done for two reasons. 
Firstly, to make sure that it contains functions having a reasonable behaviour, both, on 
the axis of symmetry and on the horizon. Secondly, such an information is needed as a 
basis for Section 5. There the resolvent of G is constructed using spheroidal functions. 
Section 6 discusses the reduced Klein-Gordon equation. The well-posedness of the initial 
value problem is shown. Further, the stability of the solutions is shown for large enough 
masses. Section 7 contains the discussion. The Appendix gives auxiliar theorems used in 
the computation of the resolvent. 




2 Conventions 

The symbols N, JR., C denote the natural numbers (including zero), all real numbers and 
all complex numbers, respectively. 



To ease understanding we follow common abuse of notation and don't differentiate be- 
tween coordinate maps and coordinates. For instance, interchangeably r will denote some 
real number greater than r + or the coordinate projection onto the open interval (r + , oo). 
The definition used will be clear from the context. In addition we assume composition of 
maps (which includes addition, multiplication etc.) always to be maximally defined. So 
for instance the addition of two maps (if at all possible) is defined on the intersection of 
the corresponding domains. 

For each \ {0}, n£N \ {0} and each non-trivial open subset M of R n the symbol 

C k (M, C) denotes the linear space of /c-times continuously differentiable complex- valued 
functions on M. Further Cq(M, C) denotes the subspace of C k (M, C) consisting of those 
functions which in addition have a compact support in Q. 



Throughout the paper Lebesgue integration theory is used in the formulation of [24 



Compare also Chapter III in [FI] and Appendix A in |50] . To improve readability we fol- 
low common usage and don't differentiate between an almost everywhere (with respect to 
the chosen measure) defined function / and the associated equivalence class (consisting of 
all almost everywhere defined functions which differ from / only on a set of measure zero). 
In this sense L 2 C (M,p), where p is some strictly positive real- valued continuous function 
on M, denotes the Hilbert space of complex-valued, square integrable (with respect to 
the measure pd n x) functions on M. The scalar product < | > on l? c (M,p) is defined by 

<f\g>--= I f*9pd n x, (8) 

J M 

for all f,g& Lq (M, p), where * denotes complex conjugation on C. It is a standard result 
of functional analysis that C*(M, C) is dense in L 2 C (M, p). 

Finally, throughout the paper standard results and nomenclature of operator theory is 
used. For this compare standard textbooks on Functional analysis, e.g.,|2B| Vol. I, [[£ 



3 Basic choices and first consequences 

As basic Hilbert space X for (|7|) we chose 

X:=L 2 c (n,g 00 <^\gt) , (9) 
where \g\ denotes the determinant of the matrix g a f,. Note that 



S°°V /r M= ( ^^ -aWflj sing (10) 

is singular at the horizon. Hence the elements of X vanish there in the mean. In the limit 
case a = this measure reduces to the standard one often used for the stability discussion 



of the Schwarzschild metric f29| , |13 |. 



The operator B is chosen as the maximal multiplication operator in X by the func- 
tion multiplying idu/dt in (|5|). Since that function is bounded and positive real- valued, 



B is a bounded linear and positive self-adjoint operator on X given by 



for every / el. The operator Aq + C is defined by 

(i4o + C7)/ := D 2 r6 f (12) 
for all / G C 2 (f2, C), where we set for every / G C 2 (Q, C) 

:= (^^-a^in 2 ^" 1 - 

(d m 2 ci 2 \ d d m 2 \ 

~fr A fr - — ~^ede sin6 d6 + ) f ' (13) 

Then Aq + C is in particular linear and (using partial integration) symmetric. Fur- 
ther again by partial integration it is easy to see that Aq + C is semibounded with the 
lower bound —a, where a := m 2 a 2 /(4M 2 r 2 _). Note that this bound approaches — oo for 
\m\ — > oo, which would suggest that the unreduced wave equation on Kerr background 
would have no stable initial value problem. Also note that it approaches for a — » 
which is the optimal bound for Schwarzschild. 



In the next step we define C := — (a + e), where e > is assumed to have the dimension 
l~ 2 . The exact value of e does not influence the results in any essential way. Finally, 
we define A as the Friedrichs extension of A . As a consequence A is a densely-defined, 
linear, selfadjoint and semibounded operator having the same lower bound as A , i.e., e. 

The objects X,A,B and C are easily seen to satisfy Assumptions 1 and 4 of p|. Ap- 
plying the results of that paper gives 

Theorem 1 (i) By 

Y := D{A 1/2 ) x X (14) 

and 

{£\0) :=< A^A^rn > + < 6b > (15) 
for all £ = £ 2 ), rj = (rji, 772) G Y there is defined a complex Hilbert space (Y, ( | )). 

(ii) The operators G and —G defined by 

G(C,v):^(-V,(A + C)C + iBr,) (16) 

for all£ G D(A) and 77 G D(A 1 ^ 2 ) are infinitesimal generators of strongly continuous 
semigroups T + : [0, 00) — > L(Y, Y) and T_ : [0, 00) — > L(Y, Y), respectively. 

(iii) For all t G [0, 00): 

|T ± (t)| ^e X p(||C|| HA-^p) ; (17) 
where \ \, \\ \\ denote the operator norm for L(Y,Y) and L(X, X) , respectively. 



(iv) For every t G R and every £ G -D(A) X /^(A 1 / 2 ) inere a uniquely determined 
differentiate map u : R — > K snca iaai 



n(t ) = e (18) 

and 

u '(t) = (19) 

for all t 6 1. i/ere ' denotes differentiation of functions assuming values in Y. 
Moreover this u is given by 

„(():={ If % (20) 
v y | T_(-t)£ fort < v 7 

/or a// t G K. 

(v) A G C is a spectral value, eigenvalue of iG if and only if 

A + C-XB-X 2 (21) 
is not bijective and not infective, respectively. 

(vi) For any X from the resolvent set of iG and any rj = (7/1,7/2) £ ^ one ^os: 

(zG-A)- 1 n = (e,^(A^ +7/0) , (22) 

w/iere 

£, := (A + C - XB - A 2 )" 1 [(£? + A) 771 - i V2 ] . (23) 

Equation (|T9|) is the interpretation of @ used in this paper. In this sense (iv) shows 
the well-posedness of the initial value problem for (H), i.e., the existence and uniqueness 
of the solution and its continuous dependence on the initial data. Moreover (p0|) gives 
a representation of the solution and (iii) gives a rough bound for its growth in time. In 
general, this bound is not strong enough to imply stability of the solutions to Part 
(v) reduces the determination of the generator's spectrum to the determination of the 
spectrum of the operator polynomial A + C — XB — A 2 , A G C jl9], |26|]. Moreover (vi) 
does the same for the resolvents. Further, || gives the following stability citeria 

Theorem 2 (i) If 

<Z\(A + C){;>+~<ti\Bti > 2 ^0. 

for all £ G D(A) with ||£|| = 1. Then the spectrum of iG is real. 

(ii) If A + C — (6/2)5 — (6 2 /4) is positive for some b G R. Tnen tne spectrum of iG is 
real and there are K ^ and £ ^ snc/i £aa£ 



/or a// 1 ^ t . 



Here | ||. | | denote the induced norm on (X, < | >) and (Y, ( | )), respectively. Note 
that the reality of the generator's spectrum would exclude the existence of exponentially 
growing mode solutions of (|5|). It seems that these criteria are not strong enough to prove 
stability of the solutions of (|5|). [] But later on (ii) will be used to conclude stability for 
the corresponding Klein-Gordon equation for cases where the mass excedes some given 
bound depending on m. Note that the positivity of A + C would imply stability via (ii). 
On first sight positivity of A + C seems unlikely because of the negative potential term 
—m 2 a 2 /A in (0). On the other hand it is well-known that the ocurrence of such a term 
can be due to the chosen representation space for A + C . In addition the domain of this 
operator is very much restricted by the condition that its elements have a compact support 
in Q. Since Q is open it follows that the support of such a function has a strictly positive 
distance from the boundary. In the theory of Schrodinger operators it is well-known from 
so called 'Poincare-Sobolev inequalities' that the kinetic energy associated with such a 
state can exceed a negative potential energy. See, e.g., |33| or for a simple example 
Vol. II, example 1 in chapter X.3. Indeed such inequalities were found, but only ones 
leading to a positive potential term with asymptotic behaviour ~ A" 13 for r — > r+, where 
^ (3 < 1. So none of them was found to be strong enough to show positivity of A + C. 
Indeed the apparent absence of better estimates lead to the impression that A + C is 
indeed negative. If this is really true it should be easy to prove using the results on the 
domain of A + C from the next section. This point has not been investigated further, 
because the negativity alone would not give any further information on the stability of 
the solutions of (El). 



4 Investigation of the domain of A+C 

In this Section the domain of A + C will be further investigated. This is done for two 
reasons. Firstly, to make sure that it contains functions having a reasonable behaviour, 
both, on the axis of symmetry and on the horizon. It turns out that this is indeed the 
case. In particular, as it should be the case, functions of the form f(r)P l m (cos6), where 
/ G C 2 (I r , C) and P l m , I = \m\, \m\ + 1, . . . , are the usual generalized Legendre polynomi- 
als are found to be in the domain of A + C . Secondly such an information is needed as a 
basis for the construction of the resolvent of G in the next Section. 



We do not give a full characterization of D(A + C) here. Instead more modestly suf- 
ficient conditions are given on functions f(r) and g(8) which guarantee that the product 
f(r)g{9) is in D(A + C). These conditions will turn out to be sufficient as a basis for the 
next Section. They are as follows: 



Theorem 3 For this denote I r := (r + , oo) and Iq := (0, tt) and define 

X r := 4(/ r ,r 4 /A) , X e := L 2 c (I dj sm6) 
and for every f G C 2 (I r , C) and g G C 2 (I e , C) 

1 



-(A/') - 



m 2 a 2 



f 



, Dig:- 



sin 6 



(sin 00') + 



m 



sin z 9 



(24) 



(25) 



1 In the following discussion the trivial cases a = 0, i.e., the case of a Schwarschild background, and 
m = corresponding to purely axial perturbations, are excluded. Of course, for these stability of the 
solutions can be concluded from Theorem 2 (ii) . 



Let be f G C 2 (I ri C) n X r and g G C 2 (I g , C) n X e such that 

D 2 r feX r and D 2 e geX e (26) 
and for m = in addition such that 

lim sin9 g' \6) = lim sax9g'{9) = . (27) 

Then f(r)g(9) G D(A + C) and 

(A + C)f(r)g(e) = D 2 r0 f(r)g(e) . (28) 
Proof: First it follows from the obvious inequalities 



r 4 (r 2 + a 2 ) 2 2 , „ 4M 2 r 4 

^Li.^^ (29) 



that L 2 c (tt, r 4 sin^/A) and X are identical as sets and that the associated norms on that 
set are equivalent. A further consequence of (p9|) along with partial integration is the fact 
that f(r)g(6) is in the domain D((A + C)*) of the adjoint (A + C)* of A + C and in 
particular that 

(A + Cyf{r)g{9) = D 2 f(r)g(6) . (30) 

Hence f{r)g(9) G D((A Q + C)*) if and only if there is a sequence h , h%,--- of elements of 
Cq(Q, C) converging to f(r)g(8) and such that for every given e > there is z/ G N such 
that for all /i, v G N satisfying fj,^ u and z/ ^ z/ : 

| < /i„ - ^|(A + C + a)(h^ -h v )>\ < e . (31) 

In the following the existence of such a sequence will be shown. Basic for this is the 
following inequality valid for all u G Cq{I t , C) and v G Cq(Io, C): 

<u(r)v(9)\(A + C + a)u(r)v(9) > ^ |w| 2 dr^ Qf sin 9v*D 2 e vd9j + 



r 4 A~Vp;! + mV/rJ>drJ (/ sintf K#)| 2 d#J < (32) 
r 4 A'V(D 2 + m 2 a 2 /r 4 +r7 2 )Mrfr ] (/ sin0v*(Z$ + l)v d9 



Here some elementary estimates have been used along with the positivity of D 2 + m 2 a 2 /r 4 h 
on C 2 (I r , C) G X r . Since A + C + a is in particular positive also the following inequality 
is valid for all Ui, u 2 G C 2 (I r , C) and t>i, t> 2 G Cq(Iq, C): 

| < wi(r)^(0) - u 2 (r)v 2 (8)\(A + C + a)[ Ul (r) Vl (9) - u 2 {r)v 2 {9)) > | 

= \\(A + C + «) 1/2 [ Wl (r) - u 2 {r)]v x {9) + (A + C + a ^ 2 u 2 (r)[ Vl (9) - v 2 {9)\ f 

< 2 < [ Wl (r) - u 2 {r)) Vl {9)\{A G + C + a) Mr) - u 2 (r)]«i(0) > + 

2 < u 2 (r)M#) - u 2 (0)]|(A) + C + a)« 2 (r)K(0) - u 2 (0)] > (33) 

where a' := m 2 a 2 /r 4 + r^ 2 . Since / is in the domain of the Friedrichs extension of D 2 on 
C 2 (I r , C) C X r there is a sequence / , /i, • • ■ of elements of Cq (J r , C) converging to /(r) 



and such that for every given e > there is Uq G N such that for all \i, v G N satisfying 
\i ^ u and z/ ^ z/ : 



r 4 A~ x (/„ - /„)• (D* + <*')(/„ - f v )dr < e. (34) 



Obviously, by an argument analogous to (|33l) this implies that the sequence 

POO 

/ r^-^+a'j/.dr.j/eN (35) 

is bounded. Moreover since g is in the domain of the Friedrichs extension of D 2 e on 
Cq(Iq, C) C Xe there is a sequence go, 9i, • • • of elements of Cq (1$, C) converging to g{9) 
and such that for every given e > there is z/o G N such that for all /x, z/ G N satisfying 
/i ^ z/ and z/ ^ z/ : 

/ sin^^-^r^l + l)^-^)^ < e. (36) 
Jo 

Here too, by an argument analogous to (|33|) this implies that the sequence 

/ sin 9 g* {D 2 e + 1) g u d9 , z/ G N (37) 
Jo 

is bounded. Finally, because of 

1 < U(r)g,(e) - Mr)g u (9)\(A + C + a)[f,{r)g,{9) - f v {r)g v {9)} > | < 

2 (J^r'A- 1 (/„ - f u )* (D 2 + a')(U - /„) dr^j (£sm9 g ;(D 2 + l)g,d9^ + 

r^- l r v (Dl + a')f v dr^ QT sintf (g„ - g„)* (D 2 e + 1) ^ - g v ) d9^j (38) 

the sequence ho, hi, . . . defined by 

K ■= Ur)g v {9) , v G N (39) 

has the required properties. □ 

In the proof there have been used facts on the Sturm-Liouville operators D 2 and Dg. 
Now, for the reader's convenience these will be given. For this define the (obviously) 
linear, symmetric and semibounded operators A r o, Ago in X r and Xq, respectively, by 

A r0 f := D 2 ,f , A eo g := D 2 g , (40) 

for every / G Co(/ r , C) and every g G Cq(Iq, C). Then one has the following 

Lemma 4 (i) A r o is essentially self-adjoint. 

(ii) Ago is essentially self-adjoint for m > 0. For m = ; the Friedrichs extension of 
Ago is given by the closure of the operator A 9F defined by A eF f := Djg for every 
g G C 2 (Ig, C)nX e satisfying flgTty together with the condition that D 2 e g G Xq. For all 
m the spectrum of the Friedrichs extension of Ag is given by {\m\(\m\ + 1), (|m| + 
l)(|m| + 2),...}. 



Proof: (i) For this define the auxiliar Sturm-Liouville operator A rQ in X r by 



A.o/:=-^(A/')' (41) 

for every / G Cq(I t , C). Obviously, A r0 is densely-defined, linear, symmetric and positive. 
Moreover since — m 2 a 2 /r 4 is bounded on I r , it follows by the Rellich-Kato theorem (see, 



e.g, Theorem X.12 in Volume II of |23|) that A r0 is essentially self-adjoint if and only if 
A r0 is essentially self-adjoint. Now, the equation (A/') = has nonvanishing constants 
as solutions. Since these are not in X r at both ends of I r , it follows that A r0 is in the 
limit point case, both, at r + and at +00. Hence A r0 is essentially self-adjoint (see, e.g., 
Finally, from this follows (i). (ii) This statement is, of course, well-known. □ 



5 Computation of the generator's resolvent 

In the following the resolvent of G will be determined for spectral parameters A which 
are non real and at the same time such that iaX is not an exceptional value.0 Note that, 
because of Theorem 1 (vi), the resolvent of G can be derived from the inverses of the 
operator polynomial A + C — XB — X 2 which are given in (ii) of the following theorem on 
a dense subset of X. 

Theorem 5 Let X be a non real element of the resolvent set of iG which moreover is such 
that iaX is not an exceptional value. For each m G Z let 

ps™(cos6>, -a 2 A 2 ) , I = |m|, \m\ + 1, \m\ + 2, • • • (42) 

be the basis [] of Xq consisting of spheroidal eigenfunctions of D 2 -\-X 2 a 2 sin 2 9 corresponding 
to the eigenvalues 

XU-a 2 X 2 ),X^ m]+1 (-a 2 X 2 ),... , (43) 
respectively.^ Finally, let be g G Co(/ r) C), m G Z and I G {|j7i|, \m\ + 1, . . . }. Then 
(i) The subset of X consisting of all finite linear combinations of elements of the form 

h{r, 6)g{r)p™{cos 6, -a 2 A 2 ) (44) 

where 

r 4 f(r 2 + a 2 ) 2 2 . 2 ^ 



h --=^[ A -a'Bm'flj , (45) 

and g, I run through the elments of Cq(I t , C) and {\m\, \m\ + 1, \m\ +2, . . . }, respec- 
tively, is dense in X . 



2 For the definition of these values see [ pp| . 

3 In the sense that the span of these functions is dense in Xq. Note that these functions are not 
orthogonal in general. Instead this sequence and the sequence consisting of its complex conjugates form 



a biorthogonal Basis of Xg. See Theorem 4 in Section 3.23 of |20 
4 For the definition of the functions ps™ see pQ| . 



(") 

{A + C-XB- A 2 ) _1 /?,(r, 6) g(r)p™ (cos 6, -a 2 A 2 ) = 
f r (r)p?(cosd,-a 2 \ 2 ) , (46) 

where f r G C 2 (I r , C) fl X r zs suc/i i/iai -D 2 / r G X r and moreover satisfies 

D 2 J r + \r(-a 2 \ 2 )(A/r*)f r = g . (47) 

i/ere /or every G C 2 (I r , C) 

Z^V := (A0')' - ^ [(ma + 2XMr) 2 + A 2 A(A + 4Mr)] . (48) 

Proof: First we notice that h is C°° on and satisfies as a consequence of ( p9|) 

r 2 /(4M 2 ) ^ /i < 1 . (49) 

Hence the maximal multiplication operator by /i in X is defined on the whole of X, 
is bijective and its inverse is given by the maximal multiplication operator Ti/h which is 
defined on the whole of X, too. by the function \ jh in X. Obviously, the subset of X 
consisting of all finite linear combinations of elements of the form 

o(rM"(cosfl,-a 2 A 2 ) (50) 

where g G Co(I r , C) and I = \m\, \m\ + 1, \m\ + 2, . . . is dense in X. Hence this is true for 
the subset of X consisting of all finite linear combinations of elements of the form 

/i(r,%(r)p™(cos#,-a 2 A 2 ) (51) 

where g G Cq(I t , C) and I = \m\, \m\ + 1, \m\ + 2, . . . , too. In the following let g be some 
element of Co(I r , C) and / be some element of {\m\, \m\ + 1, \m\ +2, . . . }. We will compute 
the element / G X satisfying 

[A + C-XB - A 2 ) f(r, 9) = h{r, 6) g(r)pf (cos 6, -a 2 A 2 ) . (52) 

We note that by Theorem 2 

(A + C-XB - A 2 ) f r {r)pf (cos 9, -a 2 X 2 ) = 

h(r,6) [{D 2 J r + Ar(-a 2 A 2 )(A/r 4 )/ r ) (r)] P T(cos6, -a 2 A 2 ) (53) 

for every f r G C 2 (I r , C) fl X r such that 

D 2 f r G X r . (54) 

In the following we construct such a f r which satisfies in particular ([|7]) Then by the 
bijectivity of A + C — XB — X 2 we conclude that 

f(r,8) = f r (r)p] n (cos6,-a 2 X 2 ) . (55) 

For this construction we need some auxiliar solutions fi, fi-, fz and of the homogeneous 
equation associated with (|47|), i.e., 



2(r-M) 

J r ' r s J r 1 



(ma + 2XMr) 2 ^ 2 / 4Mr\ s 



A 2 V A J A 



/r = , (56) 



where s := A™(— a 2 A 2 ), having special asymptotic behaviour at the singular point r = r + 
and at +00. First, by defining f r := A 1//2 / r and by introducing the new independent 
variable r* 



y/r(r + AM) + 2Mln ( (Vr + 4M + v^) 2 /(4M) ) (57) 



one gets a homegeneous first order system for f r and df r /dr* which is equivalent to 
flSlf) and which satisfies the assumptions of Theorem 4 in the appendix. From this 
theorem follows the existence of linear independent continuously different iable solutions 
(fri,df r i/dr*) and (fr2,df r2 / dr*) of the system along with continuously different iable 
functions R\ and R 2 such that 

/n(r.) = e iAr * (1 + R u (n)) , ^±(r„) = e* Ar * (iA + i? 12 (r*)) 

/ r2 (r*) = e- a ^(l+it: 21 (rO) , ^(r*) = e" iAr * HA + R 22 (n)) 
lim |#i(r,)| = lim \R 2 (n)\ = . (58) 

In the following denote by f r i, f r2 the solutions of ( p6[ ) corresponding to (f r i,df r i/dr*) 
and {f r 2,df r2 /dr*), respectively. Morover define 

f / rl for/m(A)>0 
Jrfl '~ 1 /ra for Jm(A) < . l ° yj 

Then it follows by fl58|) that 0/ rfi G C 2 (/ r ,C) n X r and D 2 r (<pf rR ) G X r for every 6 
C 2 (/ r ,R) which is identically for r < r and identically 1 for r > r l5 where r ,rx G J r 
are such that r < r%, but otherwise arbitrary. For the second step by defining g x := f r / A 
and g 2 := /' one gets a homegeneous first order system for (gi,g 2 ) which is equivalent 
to ([56]) and which satisfies the assumptions of Corollary 5 in the appendix. From this 
corollary follows the existence of linear independent continuously differentiable solutions 
(#11) #12) and (g 2 i,g 22 ) of the system along with continuously differentiable functions R 3 
and R4 such that 

gn (r) = (r - r + )- CT1 [1 + R 31 (r)} , g 12 (r) = (r - r + )^ [-i(ma + 2\Mr + ) + R 32 (r)\ 
9ii{r) = (r - r+)- (T2 [l + R Al (r)} , g 22 (r) = (r - r+j^^mo + 2AMr+) + R A2 (r)} 
lim \R 3 (n)\ = lim \R 3 (r*)\ = 0, (60) 

r»— »r+ r*— »r+ 

where 

<7i := 
cr 2 : = 



VM 2 - a 2 + z ((ma/2) + \Mr + ) /VM 2 - a 2 
VM 2 ~a 2 -i ((ma/2) + XMr + )] j^M 2 - a 2 . (61) 



In the following denote by f r3 , f rA the solutions of (|56|) corresponding to (gn,9n) and 
(921,922), respectively. Moreover define 

f ._ j Us for Im(X) > . . 

JrL • 1 / r4 for Jm(A) < . 1 > 



Then it follows by (gg) that (j)f rL G C 2 (I r ,C) H X r and D 2 r {$f rL ) G X r for every 
G C 2 (/ r .,M) which is identically 1 for r < r and identically for r > r\, where 
Tq,Ti G I r are such that r < r l7 but otherwise arbitrary. 

In the next step we notice that f r R and f r L are linear independent, because otherwise 
we would get a contradiction to the assumed bijectivity of A + C — XB — A 2 . Hence the 
Wronski determinant W of f r R and f r L 

W := A (f rL f; R - / r ' L f rR ) (63) 

is constant and different from 0. Therefore we can define 

Ur) := jf" f rL (r')g(r')dr' - / rJl (r')<Kr')dr' (64) 

for all r G J r . 

Using the foregoing results on / r£ and / r # by a simple computation It follows from the 
foregoing results on f rL and f r R and from a simple computation that / r G C 2 (I r , C) flX r , 
-D 2 / r G X r and that / r satiesfies (|47|) . Finally, from the bijectivity of A + C — XB — X 2 
we conclude (fRiP-p 

6 The case of the Klein-Gordon equation 

Compared to the wave equation considered in the previous Sections, the only change in 
this case is that the operator C has to be substituted by C := C + (ml/g 00 ), where m 
denotes the mass of the field and m^j g 00 is the maximal multiplication operator in X, 
which is defined on the whole of X as well as bounded, since this function is easily seen 
to be bounded on Q. The other objects X, A and B stay the same. Again it is easy to 
verify that X, A, B and C satisfy Assumptions 1 and 4 of 0. As a consequence one has 
theorems analogous to Theorem 1 and Theorem 2. They imply the well-posedness of the 
initial value problem, i.e., the existence and uniqueness of the solution and its continuous 
dependence on the initial data. Further, via the analogue of Theorem 2 (ii), Theorem 7 
below implies for masses satisfying (|69|), that the spectrum of the corresponding generator 
is real and that the norm of the solutions grow at most linearly in time. In particular 
there are no exponentially growing modes in these cases. 

Lemma 6 Let B' be a bounded linear and self-adjoint operator on X. Then A + B' is 
identical to the Friedrichs extension of Aq + B'. 

Proof: First, since B' is bounded linear and self-adjoint on X, it follows that 

(A + B')* = A* + B> . (65) 

Hence the domain of the Friedrichs extension (A + B')p of A + B ' is given by those 
elements / from D(Aq) for which there is a sequence /o, /i, ■ ■ ■ in D(Aq) converging to / 
and such that for every 5 > there is a corresponding such that for all v G N 

\<f»-h\{A, + B>+\\B'\){f ii -f v )> \ < 5 (66) 



if, both, /i > /i and u > u . Since fl5S[) implies 

| < - fu\A (f„ - U)> | < 5 
it follows that / is an element of D(A), too. Further, (|65|) implies 

(4, + = (A) + B')*f = Af + B'f. 



(67) 



(68) 



Hence A + B' is a linear self-adjoint (by the Rellich-Kato theorem, see, e.g, Theorem X.12 
in Volume II of |[23|| ) extension of (Aq + B')p and, finally, since (Aq + B')p is self-adjoint, 
(A + B') F = A + B'. u 



Theorem 7 Define b := ma/(Mr + ) and let be 



m ^ 



ma 



2Mr, 



2M a 2 



1 + 



+ — • 



(69) 



T/ien 

A + C + m 2 /g 00 + (b/2)B - 6 2 /4 

is positive. 

Proof: Because of the preceeding Lemma it is enough to prove that 
< f\(A + C + m 2 /g 00 + (6/2)5 - b 2 /A)f > > 



(70) 



(71) 



for all / G C 2 (Q,C). Now let / be such an element. Since its support supp(f) is a 
compact subset of Q there are r > r + and r x > r such that supp(f) C J x (0, it) C fi, 
where J := (ro,ri). In a first step one gets by partial integration, Fubini's theorem and 
Lemma 4 (ii) 



<f\(A + C)f> 
sin 9 f* 



m 2 a 2 



—A— 

dr dr A 



Id d m 2 
sin 9— + 



sin 9 89 



89 sin 2 9 



ri r* f d A d 

n -T-A- - 



m 2 a 2 



+ 



r IV 



sin^/; 



dr dr A 
1 d 



fedr 



sin # d# 

2 



sin 9 d9 
fU\m\(\m\ + r 



d in 
sin 9— + 



d9 sin 2 9 



f r d9 



dr 



ro 



m 2 a 2 



A 



fedr 



sin d9 



Further using 



< f\f/g 00 > 
<f\Bf> 

<f\f> 



/*7T r pvi 

/ / r 2 |M 2 dr 

JO Jrn 



sin 9 d9 



4mMa 



U ro 



^\fe?dr 



fdrd9 



sin 9 d9 



(r 2 + a 2 ) 2 
A 



sin9\f\ 2 drd9 



(72) 
(73) 
(74) 



we get 



< f\ (A + C + ml/g 00 + (b/2)B - b 2 /A) f > 



ri 



n h(h+i) 



m 2 a 2 r — r A 



rl r — r_ 



2„2 



m r 



2^2 



ma 



4M 2 rl 



{r 1 + 2Mr + a 2 ) f g dr 



sin Ode 



m 2 a 2 
4M 2 r+ 



fZ[\m\ + rn'{l-- 



(75) 



[(r* + 2Mr+ + 



2\ 2 2/2 

a r — r , r 



+ 2Mr + a 2 )]j / fl dr 



sin 6d6>0 



Hence it follows the positivity of A + C + mg/g 00 + (6/2)S - b 2 /4. D 



7 Discussion 

The reduced (in the angular coordinate ip) wave equation and Klein-Gordon equation 
were considered on a Kerr background and in the framework of C°-semigroup theory. 
Each equation was shown to have a well-posed initial value problem, i.e., to have a unique 
solution depending continuously on the data. Further, it was proved that the spectrum of 
the semigroup's generator coincides with the spectrum of an operator polynomial whose 
coefficients can be read off from the equation. In this way the problem of deciding sta- 
bility is reduced to a spectral problem. In addition a mathematical basis is provided for 
mode considerations. | For the wave equation it was shown that the resolvent of the 
semigroup's generator and the corresponding Green's functions can be computed using 
spheroidal functions. It is to be expected that, analogous to the case of a Schwarzschild 
background, the quasinormal frequencies of the Kerr black hole appear as resonances, i.e., 
poles of the analytic continuation of this resolvent. Finally, stability of the background 
with respect to reduced massive perturbations was proved for masses exceeding a given 
bound (see (1691)). 



It is interesting to compare the last result to earlier results of Detweiler in Q , Damour, 
Deruelle, Ruffini in H and Zouros, Eardley in ]34) . These make the existence of exponen- 



tially growing modes for the massive Klein-Gordon equation very plausible. They found 
approximate unstable modes in the superradient regime, i.e., with frequencies u satisfying 
Re(u) < ma/(2Mr + ). These modes become stable when this condition is violated. The 
approximations made in these papers lead to further restrictions. It turns out that the 
assumption of, both, these restrictions and the bound ( [69] ) derived here is incompatible 
with the assumption of superradience. Hence the stability result here does not contra- 
dict the results in these papers, but is complementary instead, f] Moreover it suggests 
that the negation of an inequality of the form of (|69"D (or some equivalent form) is the 
superradient condition. For this it should be noted that with some effort and along the 

5 To give an example for this claim, say, we would be able to show that the unstable spectrum of 
G consists of discrete eigenvalues and that the corresponding eigenstates seperate in the way assumed 
by Whiting. Then, via the results of this paper, Whiting's result |3l[] on the absence of exponentially 
growing modes would imply the stability of the solutions of the wave equation. 

6 The author is very grateful to J. L. Friedman for directing his attention to this fact. 



lines of this paper it may be possible to improve (|69|), i.e., to decrease the bound. For 
this the Poincare-Sobolev inequalities mentioned at the end of Section 3 should be helpful. 



Acknowledgements. The author is grateful to B. G. Schmidt for pointing his attention 
to the problem of defining the quasi-normal frequencies of the Kerr black hole as reso- 
nances and to J. L. Friedman, B. G. Schmidt and B. F. Whiting for valuable discussions. 

8 Appendix 

The following theorem used in Section 5 was first proved by Dunkel in || (compare also 

EBB)- 

Theorem 8 Let n G N \ {0}, a G M., I := [a, oo) and Iq := (a, oo). In addition let Aq 
be a diagonalizable complex n x n matrix and e[, . . . ,e' n be a basis of C n consisting of 
eigenvectors of Aq. Further, for each j G {1, • • • , n} let Xj be the eigenvalue correspond- 
ing to e'j and Pj be the matrix representing the projection of C n onto C.e'j with respect 
to the canonical basis ofC n . Finally, let A\ be a continuous map from I into the com- 
plex nxn matrices M(nxn,C) such that A\^ is Lebesgue integrable for each j, k G l,...,n. 

Then there is a C 1 map R : J — > M(n x n, C) with limt^^Rjkit) = for each 
j, k G 1, . . . , n and such that u : Iq — > M(n x n, C) defined by 

it 

u(t) := exp{\jt) ■ (E + R(t)) ■ Pj (76) 
j=i 

for all t G Jo (where E is the nxn unit matrix), maps into the invertible nxn matrices 
and satisfies 

u'(t) = (A + A 1 (t))-u(t) (77) 

for all t G Iq . 

This theorem has the following 

Corollary 9 Let n G N\{0}; a, to G M with a < t / H ^ N; := 1 for fi = and := /i 
for fi ytz o. In addition let Aq be a diagonalizable complex nxn matrix and e[, . . . , e' n be a 
basis ofC n consisting of eigenvectors of Aq. Further, for each j G {1, • ■ ■ , n} let Xj be the 
eigenvalue corresponding to e'j and Pj be the matrix representing the projection ofC n onto 
Ce'j with respect to the canonical basis of C n . Finally, let A\ be a continuous map from 
(a, to) into the complex nxn matrices M(n x n, C) for which there is a number c G (a, to) 
such that the restriction of A\jk to [c, to) is Lebesgue integrable for each j, feel, ...,n. 

Then there is a C 1 map R : (a, to) — > M(n x n, C) with limt^oRjk{t) = for each 
j, k G 1, . . . , n and such that u : (a, t ) — > M{n x n, C) defined by 

u(t) : = 

f E;=i(*o - t)-^ ■ (E + R(t)) ■ Pj for^ = 

I £J=iexp(Ai(fo - *)"") • (E + R(t)) ■ Pj for ^ ? ^ 



for all t G (a, to) (where E is the n x n unit matrix), maps into the invertible n x n 
matrices and satisfies 

At) = ( {to -t)^ A ° + Mt) ) ' u{t) (79) 

for each t G (a, to). 
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